Signatures of asymmetric and inelastic tunneling on the spin torque bias dependence 
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The influence of structural asymmetries (barrier height and exchange splitting), as well as inelastic 
scattering (magnons and phonons) on the bias dependence of the spin transfer torque in a magnetic 
tunnel junction is studied theoretically using the free electron model. We show that they modify 
the "conventional" bias dependence of the spin transfer torque, together with the bias dependence 
of the conductance. In particular, both structural asymmetries and bulk (inelastic) scattering add 
antisymmetric terms to the perpendicular torque (oc V and oc jelV'l ), while the interfacial inelastic 
scattering conserves the junction symmetry and only produces symmetric terms (oc n £ N). 

The analysis of spin torque and conductance measurements displays a signature revealing the origin 
(asymmetry or inelastic scattering) of the discrepancy. 
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I. INTRODUCTION 

The recent observation of current-driven magnetiza- 
tion controP in Magnetic Tunnel Junction^] (MTJs) of- 
fers promising opportunities for magnetic recording and 
memories application^^. The observed critical switch- 
ing current has now reached lO^A/cm^, which makes 
MTJs competitive candidates for Magnetic Random Ac- 
cess MemoriePI. 

However, due to the specific transport properties in 
MTJs, the characteristics of the spin transfer torque 
in these devices display significant differences with the 
current-driven torque usually observed in metallic spin- 
valveP^. Uncovering the precise form of the bias de- 
pendence of the spin torque is essential to understand 
and control the dynamical properties of the magne- 
tization. In MTJs, it has been demonstrated both 
theoreticallj^sHUl experimentally^"^ that the torque 
possesses two components, usually referred to as the in- 
plane (or Slonczewski) torque, T|| and the perpendicular 
(or out-of-plane) torque, T±. The first one is purely non- 
equilibrium and competes with the damping, whereas the 
second one arises from spin reorientation at the inter- 
faces, possesses both equilibrium (Interlayer Exchange 
Couplin^l^ and non-equilibrium components and acts 
like an effective magnetic field on the magnetization. The 
presence of this perpendicular torque results in original 
dynamical properties of the magnetizatiorP^^ti^. 

Up until now, most of the experimental efforts have 
been focused on the bias-dependence of the perpendicu- 



lar torque T±. Although this component is vanishingly 
small in metallic spin-valves, it can not be neglected 
in MTJs, due to the momentum filtering imposed by 
the tunnel barriei!^. Most of the theories, using tight- 
binding^, free-electrorPtlil or ab-initi^^ calculations, 
have addressed the bias dependence of the spin torque 
within a symmetric and purely elastic tunneling junction 
(referred to as SE tunneling). It has been shown that for 
SE tunneling at low bias voltage, the form of the spin 
torque is: 

T|| = {aiV + a2V^)Mx (MxP), (1) 
= {bo + b2V^)MxP, (2) 

where P and M are the magnetization directions of 
the pinned and free layers, respectively. These bias de- 
pendencies have been well observed in spin-diode-type 
experiments^ performed on MgO-based MTJs. The lin- 
ear bias-dependence of the in-plane torque that has been 
measured (02 ~ 0) is consistent with Ref. 112! which 
suggests that MgO-based MTJs behave like half-metallic 
junctions. 

In contrast, a number of experiments using dynam- 
ical and switching properties of the MTjJSHUl^ ^s 
well as recent theoretical investigation^^Enmi have re- 
cently questioned the validity of the "conventional" bias- 
dependencies represented by Eqs.([l)-(|2]). In particular, 
Xiao et al.^" and Wilczynski et al.HI employing the free 
electron model numerically showed that structural asym- 
metries could alter the convention bias dependence of 
the perpendicular torque, whereas Tang et al."^ predicted 
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a non-monotonic bias dependence of Tl, demonstrating 
the importance of band filhng. From the experimental 
side, Li et al!^ measured a field-like effect of the form 
cx je\V\ and interpreted their data by considering the 
electron-magnon scattering in the bulk of the ferromag- 
nets. In contrast, Sun et all^^ suggested the possibility of 
non macrospin processes or heating artefacts that would 
induce a bias-dependent effective field. Very recenlty. 
Oh et alP^ demonstrated the possibility to tune the bias 
dependence of the perpendicular torque by engineering 
the structural asymmetries of a MgO-based MTJ, con- 
sistently with theoretical simulation^^ESEI] One of the 
authors also proposed that an incomplete absorption of 
the transverse spin density within the free layer could 
lead to an asymmetric perpendicular torqu^. Finally, 
we recently studied the influence of interfacial electron- 
magnon scattering on the bias dependence of the torque^° 
and found that an additional symmetric term of the form 
oc\V\ arises. 

As seen from this brief overview, the bias dependence 
of the spin torque is far from universal, and a number of 
mechanisms has been shown to modify this dependence. 
However, the role of asymmetries has been investigated 
numerically within the tight binding modeP and the free 
electron modeP^EH and little is known concerning the 
role of inelastic scatterin^iSl20l_ this paper, we derive 
analytic solutions for T|| and T± in the case of structural 
asymmetries and (bulk and interfacial) inelastic scatter- 
ing by magnons and phonons, leading to a discrepancy 
between the actual torques and the "conventional" ones 
[Eqs. (l)-(2)]. In particular, both structural asymme- 
tries and bulk (inelastic) scattering add antisymmetric 
terms to the perpendicular torque (cx V and cx je\V\ 
), while the interfacial inelastic scattering conserves the 
junction symmetry and only produces symmetric terms 
(cx |T^|). Moreover, we suggest that a connection exists 
between the tunneling conductance and the out-of-plane 
torque which constitutes a signature of the origin of the 
discrepancy. 

This paper is organized as follows: in section II, we 
briefly discuss the form of the tunneling spin torque with 
and without spin diffusion. Section III addresses the bias- 
dependence of the spin torque and conductance in the 
presence of structural asymmetries. The influence of bulk 
and interfacial inelastic scattering is described in section 
IV and the conclusion is given in section V. 



II. SPIN CURRENT VS SPIN DENSITY 

In most of the theoretical studies on the spin trans- 
fer torque in MTJs, the torques are associated with the 
transverse spin current density at the interfaces between 
the electrodes and the tunnel barrieil^l^. This definition 
is only valid in the case of semi-infinite electrodes where 
the spin diffusion is neglected. A more correct approach 
is to relate the spin torque to the spin density rather than 
to the spin current (see Ref. dUfor a detailed discussion). 
The spin torque is then the torque exerted by the trans- 
verse spin density on the local magnetization and has the 
form: 



J 



m X Mdfi, 



(3) 



where J is the s — d exchange coupling, m is the itin- 
erant spin density, M is the local magnetization and 
is the volume of the magnetic layer. The spin density 
can be computed from the well-known spin continuity 
equatioiP^ 
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where J7s is the spin current and Tsf is the spin relax- 
ation time. In the case of a magnetic tunnel junction, 
where the resistance is dominated by the barrier, the 
spatial variation of the spin density is usually neglected 
and the torque is directly related to the spin currentPti^l 
T — — Jq^ ■ JsdO,. Therefore, in the case of a semi- 
infinite magnetic layer, the torque reduces to the inter- 
facial transverse spin current. 

However, in realistic junctions, the free layer is usually 
thin [t «2-3nm) and the torque arising at the interface 
between the barrier and the ferromagnetic electrode must 
be balanced by the torque arising at the second interface: 
T — Js{x — Q) — Js{x — t), which may introduce some 
deviations from the "conventional" bias dependence of 
the torqu^. On the other hand, in MgO-based junc- 
tions, the junction behaves like a half-metallic MTj^^J^nd 
the spin density (or transverse spin current) is strongly 
absorbed near the barrier intcrfac^^EH (a few monolay- 
ers). Therefore, the usual identification T = Js{x = Q) 
is essentially valid in MgO-MTJs if one neglects the spin 
diffusion in the electrodes. 

Nevertheless, we will show that it is possible to account 
for the spin relaxation (I/ts^ ^ in Eq. Q) in the 
bulk of the electrodes (impurities- or magnons-induced 
spin-fiip scattering) as long as this relaxation does not 
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jl^(eV) 


1.5 


1.1 


0.6 


W^(eV) 


4.67-4.81 


5 


5.04-5.35 



TABLE I: Exchange splitting and work functions for the three 
standard ferromagnetic transition metals. 



significantly modify the intcrfacial densities of states, and 
thereby the tunnehng process itself. In this case, we find 
that the resulting spin torque is a mixing between the two 
transverse components of the spin current. This issue will 
be addressed in detailed in section IV. 



III. STRUCTURAL ASYMMETRIES 



have an average exchange splitting J = {Jr + Jl)/2 with 
an asymmetry A J = Jr — Jl- To determine the influence 
of these asymmetries on the spin torque and conductance, 
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In Ref. [TSl the authors demonstrated the possibility 
to add a linear component to the bias dependence of the 
perpendicular torque by intentionally introducing struc- 
tural asymmetries in the junction. Depending on the 
asymmetry, it is possible to change the sign of the lin- 
ear component, therefore artificially tuning the form of 
the spin torque. This finding is consistent with numer- 
ical studiePiSini. Although a connection is suggested 
between the bias dependence of the conductance and the 
one of the perpendicular torque^ lo ii 18^ -j-j^jg connection 
remains unclear and analytical formulae are needed. 

In this section, we study the influence of two types 
of structural asymmetries. First, we consider the pres- 
ence of different exchange splittings in the ferromagnetic 
electrodes. The exchange splittings of Fc, Co and Ni 
have been measured experimentally near the r pointP, 
as shown in Table 1. As a consequence, varying the com- 
position of the electrodes, one can obtain different ex- 
change spHttings up to Jfl, — Ji w 0.5 eV. 

Another type of structural asymmetry is the presence 
of a different barrier height at the left and right inter- 
faces of the junction. Since the work functions of Co, Fe 
and Ni are differentt^^ (see Table |l]), the asymmetry can 
be created by using different electrode compositions, but 
also by modifying the composition of the barrier itselP^. 

We consider the junction presented in Fig. 1, where 
two ferromagnetic electrodes are separated by an insula- 
tor. The magnetizations form an angle 9 between them. 
The barrier of average height <j) = + (f>L)/2 possesses 
an asymmetry Ac/) = (pR — ipL, whereas the electrodes 



FIG. 1: Potential profile of an asymmetric Magnetic Tunnel 
Junction. The right and left parabolae represent the disper- 
sion of tunneling electrons. 

we use the same approach as Brinkman et al.^*''. Using 
the free electron model within the Keldysh formalism de- 
veloped in Ref. |2i the wave functions are determined for 
the complete structure (see Ref. IHI- The analytic forms 
of the torque and current are obtain up to the first or- 
der in exp[— 2c?Ko], where d is the barrier thickness and 
kq — ^J^jm^)]^ is the barrier wave vector for perpendic- 
ularly incident Fermi electrons (see Appendix A). The 
effective mass of the electrons within the barrier is taken 
equal to 1. Therefore, the general form of the torques 
and current is: 

T||, T_^, >/e=yy di?dk||e-2'i«(^'''ii)F(i?,k||), (5) 

where F(E, k||) is a function given explicitly in Appendix 
B, is the electron energy and k|| is the wave vector com- 
ponent in the plane of the layers. The factor £^^''''(■^'''11) 
arises from the WKB approximation and represents the 
tunneling transmission. Following the spirit of Brinkman 
et al.I^fil, we assume that the barrier is thick and high 
enough so that the energy dependence is essentially con- 
tained in the exponential factor £^^''''(■^'''11). Therefore, 
F(E,\\) « F{Ef ± eV^/2,0), and we obtain: 
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(6) 
(7) 
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at the second order in bias voltage V. The torques Tj|, 
Tl are exerted on the right layer and Gp^ap{V) is the con- 
ductance defined as Gp^apiV) = dJP^^P/dV where JP'"p is 
the charge current in the parallel and antiparallel config- 
urations, respectively. The coefficients ai.-.g^^ are given 
explicitly in Appendix C. Notice that up to the first order 
in the barrier, the angular dependence of the in-plane and 
perpendicular torques is a simple saiO. The introduction 
of asymmetries does not modify the angular dependence 
of the torque, as long as the barrier is either high enough 
or thick enough (/3 — cIkq >> 1). 

To illustrate the influence of the structural asymme- 
tries on the torques and conductance, the analytical ex- 
pressions given in Eqs. (|6|-([8| have been plotted in 
Fig. 2, together with the full numerical simulation of the 
model developped in Ref. IHl The torques and conduc- 
tance are represented in their reduced form: the in-plane 
torque is normalized to the in-plane torquance (dT^^/dV) 
in the absence of asymmetries, whereas the perpendicu- 
lar torque and conductance are normalized to their value 
at zero bias. Several points are worth noting. 

First, since the in-plane torque is already asymmetric 
against the bias voltage in SE tunneling (oi and 02 do 
not vanish in the absence of structural asymmetries), the 
"conventional" bias dependence given in Eq. ^ IS con- 
served in the presence of asymmetries and only the actual 
magnitude of ai and 02 is modified, as illustrated in Fig. 
2(a,b). Note that the small discrepancy between the nu- 
merical model (solid lines) and the analytical expressions 
(squares) can be attributed to the presence of a cubic 
term oc in the torque T| | . The change in the slope of 
the torque can be simply understood by considering the 
polarization (defined as Slonczewski's polarization - see 
Appendix B) of the electrons responsible for the in-plane 
torque T| | . 



Secondly, the terms bi and 5^'°^ are proportional to 
AJ and so that in the absence of structural asym- 
metry, the perpendicular torque and the conductance 
are quadratic in bias voltage. However, when structural 
asymmetries are present, the perpendicular torque and 
the conductance both acquire an additional linear com- 
ponent (&i and ffi '"''). This is consistent with numerical 
simulations reported earlieJ^EoiiU g^jj(j ^jjg analytical ex- 
pressions satisfactorily reproduce the numerical results, 
as shown in Fig. |2jc-f). 

An interesting feature here is the sign of the deviations. 
For A(p > and A J — 0, the junction is more conductive 
for negative bias ((/)« > 4>l), therefore a shift is observed 
in the conductance and torque towards positive voltages 
(61,5^ < - see Fig. 2(a) and Fig. 2(c)). However, in 
the case A J > and A(j) = 0, the tunneling from left to 
right is slightly more efficient (since Jr > Jl) and the 
conductance displays a shift towards negative voltages 
{g\ > - see Fig. 2(d)). In contrast, the electrons from 
the right electrodes are more polarized than the ones from 
the left electrode and the torque displays a shift towards 
positive voltages (61 < - see Fig. 2(b)). This difference 
in the signature of the structural asymmetry allows for 
the identification of the source of the linear term in the 
out-of-plane torque, as demonstrated by the study of Oh 
et al.EHl. 

These results are consistent with previous numeri- 
cal studie^^liSllil at low bias. However, the comparison 
with the tight-binding model studied in Ref. [HI presents 
some differences. The free electron model yields an 
open parabolic band dispersion whereas the tight-binding 
model produces a closed band dispersion. Therefore, the 
free electron model is only correct for low bias depen- 
dence and provides results for low hand filling. As a 
consequence, the free electron model is surprisingly well 
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FIG. 2: Bias dependence of the in-plane torque (a,b), per- 
pendicular torque (c,d) and parallel conductance (e,f) in the 
case of barrier height (a,c,e) and exchange splitting (b,d,f) 
asymmetries. The solid lines correspond to numerical cal- 
culations based on the model presented in Ref. [9] and the 
squares are calculations using Eqs. (|6|-(|8|. The parameters 
are Ep = lOeV, J = leV, (j> = 5eV, d = Inm.) 



adapted to the case of Fe. It also implies that the bias 
voltage must be smaller than the half-band width of the 
conduction electrons. As a consequence, neither band 
filling-induced sign reversal of lEC nor the oscillatory 
bias dependence^ of the perpendicular torque can be ob- 
tained within the free electron model. The above results 
are limited to reasonably small bias and low band filling 
systems. 



IV. INELASTIC SCATTERING 

In this section, we consider that the (bulk or inter- 
facial) scattering by phonons, magnons or impurities in 
the left and right electrodes are symmetric (i.e., the in- 
teractions have the same amplitude in the left and right 
electrodes). This way, the electron scattering conserves 



the symmetry of the junction (it may not be true when 
the electrodes compositions are different). Although the 
symmetry of the system is conserved, the spin torque 
does not have the same expression in the case of bulk or 
intcrfacial scattering. As mentioned in section II, in the 
case of interfacial scattering the spin torque is directly 
related to the interfacial spin current, T = J's{x = 0), 
whereas in the case of bulk scattering, the spin relax- 
ation in the electrodes cannot be neglected anymore and 
T — J's{x = 0) — J^dftm/Tsf- Therefore, although the 
symmetry of the MT J is conserved in both cases, the bias 
dependence of the spin torque will experience a different 
modification depending on whether the scattering occurs 
at the interfaces or in the bulk of the electrodes. 



A. Interfacial scattering 

We consider two types of interfacial inelastic scattering 
processes: electron-magnon and electron-phonon. The 
influence of electron-magnon scattering on TMRpZlSH and 
spin transfer torqu^^SMl has been studied within the 
Transfer Hamiltonian formalism. The current density is 
expressed in the form of a 2 x 2 spinor matrix: 

^ = 27r^ 5][pkTkp/5p(fkp)+/L(l-/fl,) 
k,p 

-/5p(Tkp)+Pkrkp/i?,(l-/L)], (10) 

where fL{R.) and Pk(p) the Fermi distribution func- 
tion and electronic density of states at the left (right) 
interfaces, and Tkp (Tpk) is the spin-dependent transfer 
matrix accounting for both elastic and inelastic tunnel- 
ing. In the spinor formalism, the charge current and spin 
current are expressed je = Tr[J] and = J's^^'' = 

Tr[a.j;{y)J], where d'x{y) are the Pauli spin matrices. 

In the case of electron scattering by interfacial 
phonons, although no spin-flip takes place, the increase 
of the conductance is expected to modify the TMPP^ 
and, correspondingly, the spin torque. In the presence of 
electron-magnon and electron-phonon interactions, the 
transfer matrix can be written: 



rpe—m 
^kp 



rpe—ph rpd 

-'kp — ^kp 




a. Si 



+ b+) I, 



(11) 



(12) 



where T^^ is the direct tunneling matrix, Q™ [Q^) is the 
phenomenological electron-magnon (electron-phonon) ef- 
ficiency, N is the number of atoms per cell, cr is the 



vector of Pauli spin matrices and Sf}^"^ are the trans- 
verse part of the magnetizations of the left and right 
electrodes. Details about the derivation of Eq. ( [TT| ) can 
be found in Ref. EQI The interaction efficiency (Q™ and 
Q^) can be related to quantum mechanical quantities, 
« J^/Ej, and Qp^" w \M^\/Ej,, where J is the ex- 
change splitting, Ep is the Fermi energy and Mq is the 
electron-phonon interactiorP^l that depends on the type 
of coupling (acoustic, optical or polar coupling). 



ously, the presence of phonons does not induce spin-flip. 
However, the direct tunneling matrix is renormalized by 

(1 + \f^^K + )) and becomes bias dependenlpl. We 
then expect that the modification of the conductance 
due to phonons alters the bias dependence of the spin 
torque. Performing the matrix products displayed in Eq. 



( 10 1 and using the definition of the spinor current stated 



above, we find: 



1. Electron- Phonon Scattering 



In the case of electron-phonon scattering, the trans- 



fer matrix Eq. ( 12 ) is diagonal in spin space and obvi- 



Je(^,q) 
T\\{E,ci) 
T±{E,q) 



G, 



e 
Go 

e 

Go 



° (1 + QtK^^) + Qt^KK)) (1 + cos0F^P^)(/i(l - fn) - fn{l - /l)). 



Pl (1 + Q^"«feq> + QJ^^^O) (/l(1 - Jr) - fR{l - II)) sine, 
(1 + Q^"(6+6q) + gP"(6q&+» (P^¥'l/l(1 - Ir) + P^^fl/fl(l - /l)) sin0. 



(13) 
(14) 
(15) 



r 



where Pl,r is the polarization at the left (right) inter- 
face and ifiR^L is a coefficient that accounts for the spin 
rotation during tunnelinJ^Hl xhe integration rules are 
described in Ref. [20l Assuming that the electron spin- 
dependent densities of state do not vary much over the 
range eV, and considering acoustic phonons (uj cx q, 
Qq oc q) with a density of states of the form pph{Lo) cx w'^, 
we obtain, at T=0 K and low bias voltage: 

G{V) - Go(l + P^P^cos0)(l + U|V^r+'), (16) 

(17) 
(18) 



II =GoP^sin^?(l-HCp/.|Fr+^)y, 
- T^o - GoP"H sm0Cph\V\ 



(ph being a coefficient that depends on the electron- 
phonon coupling, Fermi energy, Debye temperature 
Qd etc... The bias-dependence of the conductance 
(oc 1^1*^+^) is consistent with the one suggested by 
Bratkovsky when ly = 2. At larger bias, jyl""^^ is re- 
placed by ksQD and the bias dependence of the torques 
becomes linear. Note that the symmetry of the out-of- 
plane torque against the bias is conserved, whereas the 
in-plane torque acquires an antisymmetric component. 

At higher temperature and bias, more complex behav- 
ior are found, but the bias dependence of G{V) and T± 



is always an even function of y (|F|", n e N). As an 
illustration, we provide below the expressions for large 
bias at finite temperature {eV > kgT > kgQu): 



G(F) = Go(l + F^P^cos0) 1+U 



3^1 



T 



7|| =Gosin(?P^ ( 1+Cph^ (^T^ 1 \V, 



- T^o = GoP«<^L sin 05 



T ( \eV\ 



ph-. 



u+2 



Again ^ph depends on the electron-phonon coupling, 
Fermi energy, Debye temperature etc. At finite tem- 
peratures, the conductance is enhanced, due to phonon- 
assistcd tunneling and therefore, both in-plane and out- 
of-plane torques are enhanced. The temperature depen- 
dence is expected to be linear. Notice that although the 
magnitude of the torque increases with the temperature, 
its efficiency (ratio between spin torque and current den- 
sity) is not modified, since the electron-phonon interac- 
tion does not affect the spin itself, but rather the tunnel- 
ing rate. 
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2. Electron- Magnon Scattering 



In the case of electron-magnon interaction, the transfer 



matrix [Eq. (11)] possesses non-diagonal elements that 



are responsible for spin-flip scattering. We then expect a 
much more complex influence on the torque. Assuming 
a magnon density of states of the form Pm{^) = i 
symmetric electrodes {P^ — — P, (pL — iiR — 4>) 
and T=0 K, we find: 

G(y) cx (l-p2cos6l)|l/r+\ (22) 
Til - I||o oc sin6l[P(l + P) - (1 - P)(l + Pcos6l)]F'(^S) 
Tl - Tlo cx P0sin6'(l - coBe)\V\''+'^ . (24) 

The detail of these expressions can be found in Ref. 
1201 Interestingly the perpendicular torque and the con- 
ductance both acquire a component that is symmetric 
against the bias. Furthermore, since the electron-magnon 
interaction mixes the majority and minority channels, the 
angular dependence is also affected, contrary to the case 
of electron-phonon coupling. 

The finite temperature situation has been studied in 
Ref. 20 and gives rise to a non-linear dependence as a 
function of both voltage and temperature. Actually, com- 
peting mechanisms take place when both magnon emis- 
sion and absorption are accounted for. Let us consider 
the torque exerted on the right electrode magnetization. 
Magnon emission (absorption) occuring at the left in- 
terface increases (reduces) the effective spin-polarization 
of the incoming electrons, therefore enhancing (lowering) 
the spin torque exerted on the right electrode. Symmetri- 
cally, electron-magnon interactions occuring at the right 
interface also affects the effective polarization of electrons 
coming from the right reservoir. Finally, we must stress 
out that the detailed temperature and bias dependencies 
presented here are strongly conditioned by the electrons, 
phonons and magnon band structures. 

B. Bulk scattering 

In contrast with interfacial scattering, in the case 
of bulk scattering (by impurities or magnons) the spin 
torque is no more described by the purely interfacial spin 
current since spin relaxation can not be neglected in the 
bulk of the layers. Therefore, the spin torque reads: 



J 



T = / X MdVL = / [-VJs - —]dVl (25) 

The presence of a finite spin relaxation time {1/Tsf ^ 0) 
induces a coupling between the two components of the 



spin torque, so that the perpendicular torque now in- 
volves a contribution of both in-plane and perpendicular 
interfacial spin current densities. In a MTJ, the inter- 
facial densities of state are usually only affected by the 
first few monolayers away from the interface. Therefore, 
since the spin-diffusion length is on the order of 5-15nm, 
we can assume that the tunneling process is almost not 
affected by the presence of spin-flip scattering. Then, 
the interfacial spin current can be identified to the spin 
torque without spin-flip: J7s(x = 0) = Tq. As a conse- 
quence, the actual spin torque has the form: 



^11 



T\\o ~ 



Tj 
T.J 



(26) 
(27) 



where tj = h/J. In the case of a symmetric magnetic 
tunnel junction in the absence of interfacial inelastic scat- 
tering, T||o and T±o are given by Eqs. ([I]) -([2]). For low 
bias voltage, when the spin-flip is dominated by EUiott- 
Yafet spin scattering, r^f is bias-independent (but tem- 
perature dependent) and the perpendicular torque gains 
a linear component aiVrj/rsf. 

At large bias, or non-zero temperature, the spin-flip 
scattering is dominated by the electron-magnon interac- 
tion. As showed by Li et al^, the spin-flip relaxation 
time due to electron-magnon interaction is inversely pro- 
portional to \V\. Consequently, the presence of bulk 
magnons results in an additional component in the per- 
pendicular torque of the form oc je|^|- This is in sharp 
contrast with the case of interfacial magnons, where the 
additional component is simply \V\. 

V. CONCLUSION 

In summary, we studied the influence of structural 
asymmetries and inelastic tunneling on the bias depen- 
dence of the spin transfer torque in MTJs, using either 
the free electron model or the Transfer Hamiltonian for- 
malism. Our results are summarized below: 

1. Structural asymmetries: the perpendicular torque 
and the conductance acquire a antisymmetric linear 
component of the form oc V, while the bias depen- 
dence of the in-plane torque is still described by 
Eq. ([T]). The obtained formulae provide consistent 
results in the low-bias region and at low band fill- 
ing with the numerical results of the tight binding 
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modeP and are in good agreement with the numer- 
ical results of the free electron modetiSliiJ. Conse- 
quently, they can serve as a guideline to design the 
spin torque bias dependence, as demonstrated by 
Oh et alP. 

2. Inelastic interfacial scattering: the symmetry of 
the MTJ is conserved and the perpendicular torque 
and conductance acquire a symmetric linear com- 
ponent of the form cx n E N. The presence of 
magnons or phonons interactions is usually revealed 
through peaks in the conductance derivative. The 
influence of the temperature has been briefly dis- 
cussed. 

3. Bulk spin-flip scattering: the spin torque is no 
more equal to the net transfer of angular momen- 
tum. The relaxation of the spin accumulation in- 
duces a mixing between the two components of the 
torque, giving rise to an antisymmetric component 
of the form V and je\V\ in the case of impurity- 
and magnon-induced spin scattering, respectively. 
Since the resistance is dominated by the barrier, the 
contribution of bulk scattering is usually negligible 
on the conductance. 

Finally, we suggest that a link exists between the sig- 
nature of asymmetry and inelastic scattering in the per- 
pendicular torque and conductance. Since both are sym- 
metric against the bias in a symmetric MTJ, the intro- 
duction of structural asymmetries or inelasticity affects 
both quantities, but in different ways. The careful anal- 
ysis of the perpendicular torque together with the con- 
ductance should give important clues on the origin of the 
additional linear terms, as suggested in Ref. US] in the 
case of structural asymmetries. Note however, that the 
conductance remains unaffected by bulk scattering and 
therefore, the influence of bulk magnons cannot be an- 
alyzed by comparing the perpendicular torque and the 
conductance. 



Appendix A: Wave functions in the large barrier 
approximation 

We use the free electron model within the Keldysh for- 
malism as described in Ref. [9l The electron wave vectors 
for majority and minority spin in the left and right elec- 
trodes and in the barrier are then: 



fci^2 = y|?(^^~£^ll±^.-f ) (Al) 
fc3,4 = ^E-E^^±Jn + f^ (A2) 

(0L + ^+i?F-i? + i?||-^(eF-4<^)) 

The indices 1,3 (2,4) refer to the majority (minority) spin. 
The wave function of an electron of injected from the i- 
th electrode with an initial spin a is represented in the 
vector form = (^V' ^lo"-'' ^^^^ wave functions for 
the electrons from the left and right electrode at the in- 
terfaces are therf^ 





V2h 

kl + ZKl 


(A4) 




den 


sine (A5) 




k2 + ini 


(A6) 


n 


den 


sine (A7) 




= AiEn (^2 + 

den 


iKi){k4 + iK^) cos ^A8) 




= AiE„ ^/2^3«:iK2 ^ + j^^) sin ^9) 
den 2 




%/2fc4KlK2., , 

= AiEn -J (fci + 

den 


iKi){k^ + ii^2) cos^lO) 




= 4zS„ (k2 
den 


Q 

+ iKi){k-i + iK.2) si(i«^l) 
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with den — 2E'^{ki + iKi){k2 + iKi)(fc3 -|- iK2){k4 + iK2), 
En — exp — (iy^^r|i ndx is the exponential factor and 
k(x = 0) = ki^k(x = d) = K2. The above equations 
together with the integration rules mentioned in section 
II. A. are sufficient to describe the transport properties of 
the junction. 



Appendix B: Currents and Torques in the large 
barrier approximation 

By definition, the charge and spin currents are defined 

as 

J» = ^ y J dEdkii ((a, ® V)l/l + (ct, ® V)fl/H) (Bl) 

where i = 0,x,y, z and <t = {ux,Oy, Oz) are the spin Paul 
matrices, ctq is the identity and {■■.)l,b. denotes quantum 
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mechanical averaging, involving the rightward and left- 
ward spin-dependent wave functions defined in Ref. IHl 
and ffj are the Fermi distribution functions of the left 
and right reservoirs. Expanding these wave functions up 
to the lowest order in the barrier height, the charge and 
spin currents for a majority electron issued from the left 
reservoir are: 



dEdki 



iKiK2(fc3 -t- kijinf + fc|)(K2 -I- k^ki) 



ini + kl)ini + ki)ini + ki)i4 + ki) 



[1 + Pl cos 6] h 



dEdki 



4fciKiK2(fc3 - k4){Kl + kDinl - k^ki) 



dEdk 



{n\ + kl){Kl + kl){Kl + kl){4 
'ikiKiK\{k1 — k\){K\ 



{n\+kl){Kl + kl){nl 
I 



k1){n\ 



kl) 



Jl sin 9 



fh sin e 



(B2) 
(B3) 
(B4) 



and Pl 



(fci-fc2)(Ki-fcifc2) 



is Slonczewski's 



(fcl-|-fe2)(K?-|-fcl/C2) 

polarizatiorii^.The contribution for a minority electron 
is obtained by performing the following replacements: 
^1,3 ^ ^2,4 and 9 — >■ —9. Similarly, the contribution 
of electrons issued from the right reservoir is obtained 
by performing the following replacements: Ki -f-^- 
(1,2) o (3,4) and — > Jr. The final expressions are 
then: 



T _ tT 1 _ tT _ 7^ 

■^e — "JeL ' -^eL '^eR eR 



(B5) 



Til =T 



Ml 



^\\R 



^\\R 



(B6) 
(B7) 



Appendix C: Analytical expressions for Spin 
Torques and Conductance 



After some algebra using Eqs. (|6|-([8]), 
following results: 



we obtain the 



T\\o 



^0 



{k'i 



fc?)(KQ 



2m 2n^/3 (k§ + fc2)2 («•,§ + fc2)2 



nl {k^^-kl){k^ 



-2P 



k'^kj) 



-2/3 



2m 7r2^2 



klY 



2e^ 
h 7r2/3 



(Kq -I- fc2)2(K2 , ^ 
-k^f{k^k^ + Klf±{k^ 



fc2)2(, 



k]Y 



^2\2 



-2/3 



(CI) 
(C2) 
(C3) 



and 



Kl[2{l3-l)k^kl + Kl{e^ + kl)] AcP 

2k^k^{{l3 - I)fc2fc2 _ (;3 _ 3)^4 _^ (fc2 + ^2)^2) ^ 
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02 



P 

9i 



ap 

9i 



92 



8fc^fc^((/3 - l)^kl - (/3 - 3)4 + + kl)4) 



4fc3fc3((/3 - l)/c|fc| - (/3 - 3)k^ + (fc| + kl)Kl) 

(4 



kinw- 



l)fcl 



AJ 



/3 



Sk^kfii^- 
2 k^kiik^ki 



l)fc|fc| 
3«;2)(/c^ - 



- 3)k;q 4- (fc^ 1 n-^ 

(2/c^fc;-K2)(^3 



8 

k^k^iP + 3) 



(2/3 - 9)Kg 



fci)('«o - k^ki) 



6{k^ki 
5 

~ 12, 



klfAJ 



(«2 + fc2)2. 



12 



3\ A(j) 



fc2(K2 + fc2)2 
f k2)(/c2 + ^2) ^ ;.2(3;.2 _^ ^2)(^2 ^ ^2)] 



AJ 



16fc4fc4(K2 



fc2)(«:2. 



J 



^A0 
12 6 



(C4) 



(C5) 

(C6) 
(C7) 
(C8) 

(C9) 
(CIO) 



^(1). When 



where k^^^ — y^(£'i?± J) and kq - 
the barrier becomes thinner, a corrective multiphcation 
factor of the form (1 + ^ + 4^2 ) should be inserted into 



Eqs. (C1|-(C3). Note that Tj_q and Go are similar to 



previous derivations using a free electron modeP^Ml. The 
above relations are limited to low bias voltage in low band 
filling systems. Using more realistic densities of states, 
these relations may by modified. 
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